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f**^ ' The recent type la supernova data suggest that the Universe is accelerating now and 

' decelerated in recent past. This may provide the evidence that the standard Friedmann 

, equation needs to be modified. We analyze in detail a new model in the context of mod- 

ified Friedmann equation using the supernova data published by the High-z Supernova 
Search Team and the Supernova Cosmology Project. The new model explains recent 
acceleration and past deceleration. Furthermore, the acceleration of the expansion of 
^ , the Universe is almost zero in the future. 
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1. Introduction 



The type la supernova (SN la) observations support that the expansion of the 
Universe is speeding up rather than slowing down ^'^-S. A spatially flat universe is 
favored by the measurements of the anisotropy of the cosmic microwave background 
(CMB) '^'^ The observation of type la supernova SN 1997ff at z ~ 1.7 also provides 
the evidence that the Universe is in the acceleration phase and was in the decel- 
eration phase in recent past The transition from the deceleration phase to the 
acceleration phase happens around the redshift zt ^ 0.4 The simplest way of 
explaining the acceleration phenomenon is through the cosmological constant with 
the equation of state parameter lua — ~1- One easily generalizes the cosmological 
constant model to dynamical cosmological constant models, such as the dark en- 
ergy model with negative constant equation of state parameter — 1 < loq < —1/3 
and the holographic dark energy models In general, a scalar field Q called 
quintessence field that slowly evolves down its potential V{Q) takes the role of a 
dynamical cosmological constant 11.12,13,14^ r^^ie energy density of the quintessence 
field must remain very small compared with radiation and matter at early epoches 
and evolves in a way that it started to dominate the Universe around the redshift 
0.4. If we remove the null energy condition restriction uq > — 1 to allow superneg- 
ative ujQ < — 1, then the phantom energy models arise More exotic equation of 



1 
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state is also possible, such as the Chaplygin gas model with the equation of state 

p = —A/p and the generalized Chaplygin gas models with the equation of state 
p = —Ajp'^ . There are fields other than the quintessence field, such as tachyon 
field, as dark energy 

Although there arc many dark energy models proposed in the literature, the 
nature of dark energy is still unknown. Therefore it is also possible that the observa- 
tions show a sign of the breakdown of the standard cosmology. In brane cosmology, 
our observable universe is a 3-brane which is located at the symmetry fixed plane 
of a five dimensional Ads (Anti-de Sitter) spacctime with the extra dimension to 
be an orbifold SxjZi. Binetruy, Defi^ayet and Langlois showed that the Friedmann 
equation takes the unconventional form B."^ ~ p + in brane cosmology In 
addition, Chung and Freese argued that almost any relationship between H and p 
was possible with the help of extra dimension Along this line of reasoning, the 
usual Friedmann equation = 8ttGp/3 is modified to a general form = g{p) 
20,21,22,23,24,25,26 ^jiggg modcls, HO exotic matter form is needed and the Uni- 
verse is composed of the ordinary matter and radiation only. The function g{x) 
satisfies the following conditions: (1) 5(2:0) = 1; (2) g{x) « x when z ^ 1; (3) 
g{xQ) > 3xog'{xo)/2, where x = fimo(l + z)^ + ^^ro(l + z)"^, Xq = OmO + ^rO ~ ^mO 
and firo = 8.40 x 10~^ ^. In Cardassian model, g{x) = x + Bx^. In generalized 
Cardassian model, g{x) = x[\ + Sa;"*^"^^^]^/" For Chaplygin gas like model, 
g{x) = X + {A' + B'x^y^^ In Ref. 21, Friedmann equation is modified to be 
jj2 _|_ jja ^ p based on the idea of embedding a three brane into a five dimen- 
sional flat spacetime. The modified model can be mapped into dark energy model. 
For example, g{x) = x for standard cosmology; g{x) = .t + A for ACDM model; 
g{x) = x + iJx^+'^'S for dark energy models with constant equation of state param- 
eter wq. In this paper, we discuss the new model proposed in Ref. 26 in detail. We 
first derive the model with the help of extra dimension, then we use different SN 
la data sets to fit the model. 

2. Model Setup 

Follow Chung and Freese the new model motivated by extra dimensions can be 
derived from the action 



where M is the Ricci scalar in five dimensions, Q is the five dimensional metric 
determinant, K5 = ^/rripi'^^ is the five dimensional Newton's constant. Sorb is the 
orbifold or the bulk action, and ^boundary represents the boundary or the orbifold 
fixed plane action. In this setup, our universe is a 3-brane residing at the boundary 
of the fifth orbifold SxjZi. Varying the above action with respect to the metric, we 
get the five dimensional Einstein equation 




(1) 



GmN = ^MN 



(2) 
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where Tmn is the energy momentum tensor of the bulk and Imn is the energy 
momentum tensor of the boundary. 

To get cosmological solutions, we take the general five dimensional metric for 
an isotropic and homogeneous flat brane universe embedded in five dimensions as 

ds^ = -N'^{t,y)df +a'^{t,y)Sijdx'dx^ +B'^{t,y)dy'^. 

The components of the five dimensional Einstein tensor are 



(3) 
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(7) 



where a = da/dt, a = (Pa/dt^ and a' = da/dy. Assume that the matter in the 

visible brane located at the fixed plane y = takes the form of the perfect fluid, 
then we get = 5{y)dia.g{—p,p,p,p, 0)/B. Evaluate Eq. (2) for the Gqo component 
Eq. (4) and the Gu component Eq. (7) at the boundary y = 0, we get the junction 
conditions 



aOBO 



[N'] 



4 



-(2p + 3p), 



(8) 



(9) 



NOBO 3 

where superscript means that the variable takes value at y = and the junction 
[a'] = a'(t,0+) — a'{t,0-). The junction of the G05 component Eq. (6) gives the 
usual energy conservation equation 

d" 

P + 3— (p+p) -0. 



The average of the G55 component Eq. (5) at y = together with Eqs. (2), (8) and 

(9) gives 



(a°) 



0\2 



(iV0)2 1(^0)2 a^NO a'^l 3(B0)2""" 35 
In terms of the cosmic time dr = N^dt, the above equation can be rewritten as 



(10) 
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Now if we fix the gauge = 1 and assume that T55 is a constant, then we get the 
four dimensional Friedmann equation 

where C is an integration constant. If our 3-branc is embedded in a five; dinicmsional 
Ads spacetime, then T55 = and p = T + pb, here — A^ is the five dimensional 
cosmological constant, pb is the matter density in the brane and T is the brane 
tension. If for some reasons, we have fine-tuned = QAb/k^ and C = 0, then from 
Eq. (12), we get the familiar brane cosmological equation 

^ -36^"+ 18^^^- 
If the bulk energy momentum tensor takes the perfect fluid form = 
dia,g{—pB,PB,PB,PB, —Pb), then we can get the solutions to Eq. (2) as follows: 

B{t,y) = l, (13) 



ATI, \ (14) 

a^{t, y) = A{t) cosh{py) + B{t) &mh{p\y\) + C(t), (15) 



where pb = 2k|/9b/3 in the case that pB < and C{t) = "ia^it)/ {k\pb), or 

a\t, y) = A{t) cos{py) + B{t) sm{p\y\) + C{t), (16) 



where pb = \/2i^Pb/^ in the case that pB > 0, or 

aHt,y) = a^{t)y^ + D{t)\y\ + E{t), (17) 

in the case that pB = 0. Furthermore, we can fix the temporal gauge such that N'^ = 
1. The integration constants A{t), B{t), D{t), E{t) and a{t) can be determined from 
the junction condition Eqs. (8) and (9). 

If the bulk energy momentum tensor takes other forms, Chung and Freese 
argued that almost any relationship between p and H was possible in Ref. 19. To 
see this, we take the matter source in our brane Universe to be dust, i.e., p = 
Pb oc {aP)~^ and p = pb = 0. Then the junction condition Eqs. (8) and (9) can be 
rewritten as 

2M __JiV^__2«| 
aOfiO - Tvo^o - 3 Pb- 

If we need a relation pb = uH'^, then we shall look for the solution to the following 
equation 

The solutions to Eqs. (18) and (19) are 

N{t,y) = B{t,y)=eMm\y\), (20) 
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a{t, y) = '' exp(-/3(t)[F(|y|) - F{y = 0)]/2), (21) 

where (3{t) = K.'^u{q/3t)'^ /3, F{y) is any smooth function satisfying F'{y = 0) = 1 
and to is an integration constant. Substitute the above solutions to Eq. (2), we can 
get the bulk energy momentum tensor. 



2.1. New Model 

In ref. 26, we proposed a new Friedman equation 

F2 = ap6[l + exp(-7p,)]". (22) 

In this section, we show how to derive the model. Now wc arc looking for the 
solutions to the full five dimensional Einstein equation (2) satisfying the boundary 
condition Eq. (18) and the following equation 

The general solutions which satisfy Eqs. (18) and (23) are 

N{t,y) ^ B{t,y) = cMKlPb{t)\v\/i), (24) 

a(i, y) = a°it) cxp{- Klph{t)[F{\y\) - F{y = 0)]/6), (25) 

where F{y) is any smooth function satisfying F'{y = 0) = 1, pb{t) oc (a°)~^ and 
aP{t) are the solutions to Eq. (22). Substitute the above solutions into Eq. (2), 
we can get the five dimensional energy momentum tensor which provides the new 
model. 

In general, the modified Friedmann equations for a spatially flat universe are 

= Hlg{x), (26) 

p + 37?(p + p) = 0, (28) 

where x = SttGp/SHq, 1 + z = ao/a is the redshift parameter, a subscript means 
the value of the variable at present, g{x) is a general function of x and from now on 
g'{x) = dg{x)/dx. Let fJmO = 87rGpmo/3-ffo' then xq = fimo and x = Oi„o(l + z)^ 
during matter dominated era. For our new model, we have 

gix) = .t(1 + e""^)" = a;[l + {x^^^" - 1)^/^°]", (29) 

where a = — ln(f2~Q^" — l)/riniO > 0. 

Dark energy models have been extensively studied in the literature. Most re- 
searchers studied dark energy models by constraining the dark energy equation of 



(23) 
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state parameter cuq. To compare the modified model with the dark energy model, 
we make the following identification 

xg'jx) - g{x) 

WQ = ;r- . (30) 

g[x) — X 

The transition from deceleration to acceleration happens when the deceleration 
parameter q = —a/aH'^ = 0. From equations (26) and (27), the transition redshift 
Zt is given by 

g[n^o{l + z^f] = ^0^o(l + ^T)'5'[^^mo(l + ^t)']. (31) 

3. Supernova fitting result 

The current age of the Universe is 

"""'"^""Vo ''^io (l + %V2[0.o(l + .)3] 

The luminosity distance c^l is defined as 

rfL(^) = ao(l + z) r ^ = ^^^^ r g-^/'[ClM^ + uf]du. (33) 

Ji a(t ) Ho Jo 

The apparent magnitude redshift relation is 

m{z) = M + 5 logio dL{z) + 25 = M + 5 log^o T^i^iz), (34) 

where 'Di^{z) = Hodi^{z) is the "Hubble-constant-free" luminosity distance, M is 
the absolute peak magnitude and = M — 5 log^o -ffo + 25. The parameters in our 
model are determined by minimizing 

[mohs{zi)- m{zi)Y 
a. 



where oi is the total uncertainty in the observation. The ^^-minimization procedure 
is based on MINUIT code The observational data used are: the 20 radio galaxy 
and 78 supernova data by Daly and Djorgovski the 172 supernova data with 
z > 0.01 and Ay < 0.5 mag by Tonry et al. and the 22 supernova data with 
Ay < 0.5 mag by Barris et al. ^, and the latest 157 supernova data compiled by Riess 
et al. ^ combined with the CMB shift parameter TZ = i}]^^ Hodi,{zis) / {c{l + zis)) = 
1.716 ± 0.062 here z^ = 1089 ± 1 ^ For the SN la data compiled by Tonry et 
al. and Barris et al, a 500 km/s uncertainty in quadrature to the redshift errors is 
added to account for the velocity uncertainty and is calculated by marginalizing 
over Hq. The marginalization is done by integrating e~^ over all possible values 
of Hq. Since Hq is appeared linearly in x^, the integration process is equivalent to 
minimize over Ho- 

From equation (29), it is easy to see that at high redshift z, x is very large 
and g{x) ~ x. Therefore the standard model is recovered at early times. In the 
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future, exp{—ax) 1 and ^(a;) « 2"a;, so the Universe will again evolve as t^^^. 
Substituting the new model (29) into equation (27), during the matter domination 
we get 

= -lx{l + e-"^)" + ^nax^il + e-«^)"-ie-«^ (36) 

As a; 1, we have a/{aHQ) —x/2. As a; ^ 0, d/{aHQ) — > 0. In this model, 
equation (31) becomes 

1 + (O-J/" - l)(i+-)^ + 3nln(0-J/" - 1)(0-J/" - l)(i+-)^(l + .t)^ = (37) 

If we compare this model with the dark energy model, from Eq. (30), we get 

^ n(l-OL/o")ln(0-;/"-l) _ 

f^mo(l — f^mo) 

In the fit, the range of parameter space for n is (0, 45] and the range of param- 
eter space for fimo is (0, 1]- The best fit results are summarized in table 1. In table 
1, we also list the transition redshift, the state of equation parameter ujqo of the 
equivalent dark energy model and the age of universe using the best fit parameters. 
For the last row in table 1. we fit the model to both the 157 gold sample supernova 
data and the CMB shift parameter. In this model, the Universe is older than that 
by the standard cosmology. The contour plots to different data sample are shown 
in Figs. 1-5. 
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In Fig. 4, we plot the dimensionless acceleration a/aH^ against the redshift. 
It is clear that the model have deceleration phase in recent past and acceleration 
phase at present. In Fig. 5, the evolution of wq of the equivalent dark energy 
model is shown for the best fit parameters to the combined supernova and CMB 
data. Around z r-^ 2, g{x) x and the Universe became indistinguishable from the 
standard matter dominated universe. 

As shown in Fig. 4, the new model gives past deceleration, current acceleration 
and future deceleration or near zero acceleration. From theoretical point of view, 
future deceleration is favored because it avoids the horizon problem. In Ref. 10, the 
authors show that eternal accelerating dark energy models prevent us from ever 
measuring inflationary perturbations which originated before the ones currently 
observable due to the growth of the comoving Hubble scale in the future. The 
Universe in this new model is older than that in standard cosmology as shown 
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Q. „ and n Confidence Contour 

mO 




Fig. 1. The OmO and n contour plot for the 20 radio galaxy and 78 supernova data by Daly and 
Djorgovski 




Fig. 2. The OmO and n contour plot for the 194 combined supernova data by Tonry et al. and 
Bajrris al. 

in table 1. Gong and Daly and Djorgovski found the model independent result 
zt ^ 0.4 ^. The result is consistent with ours as seen from table 1. The somewhat 
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Fit to Daly & Djorgovski data 
Fit to Tonry & Barris et al. data 
Fit to Riess et al. & WMAP Data 




Fig. 3. The Icr QmO and n contour plots 




Redshift z 



Fig. 4. d/aH'^ versus tlie redsliift z for tiie best fit parameters listed in Table I. 



model independent study by Alam, Sahni, Saini and Starobinsky found that wq 
evolves rapidly from zero at high z to a strongly negative value at present This 
result is also consistent with ours as shown in table 1 and Fig. 5. In conclusion, we 
proposed a new model which is consistent with current supernova observations and 
avoids the future horizon problem. It is difficult to discriminate this model from 
dark energy models with current observations. Future supernova data from the 
dedicated satellite telescope (such as Supernova/ Acceleration Probe) is necessary 
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Redshift z 



Fig. 5. U!Q versus the redshift z for the best fit pajrameters HmO = 0.52 and n = 2.8. 
to differentiate these models. 
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